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Abstract 

Variational calculation of the ground state energy and its properties using the second-order 
reduced density matrix (2-RDM) is a promising approach for quantum chemistry. A major obstacle 
with this approach is that the iV-representability conditions are too difficult in general. Therefore, 
we usually employ some approximations such as the P, Q, G, Tl and T2' conditions, for realistic 
calculations. The results of using these approximations and conditions in 2-RDM are comparable 
to those of CCSD(T). However, these conditions do not incorporate an important property; size- 
consistency. Size-consistency requires that energies E(A), E(B) and E(A - ■ ■ B) for two infinitely 
separated systems A, B, and their respective combined system A - ■ ■ B, to satisfy E(A ■ ■ ■ B) = 
E(A) + E(B). In this study, we show that the size-consistency can be satisfied if 2-RDM satisfies 
the following conditions: (i) 2-RDM is unitary invariant diagonal iV-representable; (ii) 2-RDM 
corresponding to each subsystem is the eigenstate of the number of corresponding electrons; and 
(iii) 2-RDM satisfies at least one of the P, Q, G, Tl and T2' conditions. 
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I. INTRODUCTION 



We are interested in direct determination of the second-order reduced density matrix (2- 
RDM) as an alternative to using the wavefunction, because 2- RDM contains all the relevant 
information for the ground state energy and physical properties of the iV-particle (electronic) 
system in a compact form Qj. In 2-RDM formulation, the ground state energy is obtained 
by minimizing the Hamiltonian, which is a linear functional of 2-RDM. We call such an 
approach as "the RDM method" . 

The iV-representability conditions in 2-RDM Q are the key to the RDM method and 
essential for electron correlation. However, complete establishment of the 2-RDM seems 
impractical to date (3, ^]. As such,we usually employ and investigate some of the necessary 
conditions such as the P, Q j^], G [3], Tl, and T2 5, 6] conditions. All these linearly positive 
semidefinite conditions involve the 2-RDM; Erdahl and Jin have previously extended these 
representability conditions that involve the higher-order RDMs {j]. 

The first calculation employing the P, Q and G conditions as approximate N- 
representability conditions was done by Garrod et al, and they have reproduced very good 
ground state energy for the Be atom ja, Q]. However, when Mihailovic et al. applied these 
conditions to nuclear structure calculation [l(J with disappointing results, researchers there- 
after seemed to lose their motivation. 

About 25 years later, in 2001, Nakata et al. formulated the problem in a primal semidef- 
inite program, and performed direct variational calculation of the 2-RDM employing P, Q 
and G conditions using a well-established semidefinite programming solver called SemiDef- 



inite Program ming A 
and molecules 
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gorithm (SDPA) 11| and applied the method to many small atoms 
13]. Their results with the P, Q, and G conditions were very encour- 
aging; yielding 100 to 130% correlation energy, and they even reproduced the dissociation 
curve of the nitrogen dimer in good agreement with fullCI. When Zhao et al. subsequently 
reformulated and applied the Tl and T2 conditions to small molecules their results 
are comparable to CCSD(T) at equilibrium geometries. The calculated correlation energies 
range typically from 100 to 101%, and these results are evidently better than using P, Q 



and G conditions 
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13] and Zhao et al. 



15] . Mazziotti immediately confirmed the findings of Nakata et al. 
5] , and applied these conditions to larger systems by implementing 



a semidefinite programming solver 
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18] . An attempt by Mazziotti et al. in massive ap- 



2 



plication to verify these iV-representability conditions is in fact important. Later, Braams 
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20]. 



et al. and Mazziotti found T2' condition, which slightly improves T2 condition 

Therefore, necessary conditions such as P, Q, G, Tl and T2' seemed very promising; 
however, in 2009, Van Aggelen et al. reported a drastic failure of the dissociation limits 
of some molecules on dissociation 21|. As they observed fractionally charged atoms, they 



become aware of P, Q and G conditions lacked size-consistency. For example, at the disso- 
ciation limit of CN~, the Mulliken population of C and N~ were 6.60 and 7.40, respectively. 
When Nakata et al. examined size-consistency of 32 non-interacting methane and nitrogen 
molecules under P, Q and G conditions, non-size-consistent contributions were 3 x 10 -4 



and 3 x 10 3 atomic unit 



221 ] . Verstichel et al. proposed a method to fix these patholog- 



ical behaviors by imposing fractional iV-representability conditions to the subsystems with 
numerical verifications, an d ap plied them to the dissociation limit of BeB + , N2 and isoelec- 



tronic molecule series 
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241 ]. Although useful in some cases, this method would not fix 



non-size-consistency problems for the all cases. Moreover, calculation of various electron 
numbers for the subsystems is required, and the choice of subsystems is unclear, especially 
when a chemical reaction occurs 22]. 

In this study, we show condition adequate for establishing size-consistency; an approx- 
imately iV-representable variational space is size-consistent when 2-RDM satisfies the uni- 
tary invariant diagonal iV-representability, and positive semidefinite type iV-representability 
conditions, assuming that the calculated 2-RDM for each subsystem is an eigenstate of the 
number of electrons. The diagonal iV-representability conditions of 2-RDM have been ex- 
tensively studied by Yoseloff, Kuhn, Davidson, McRae and Pistol 25M27I]. Useful features of 
these conditions are that they can be written only by linear inequalities of the diagonal ele- 
ments of 2-RDM, and if a given 2-RDM is diagonal ^-represent able then there exists certain 
wavefunction that reduces to the diagonal elements of 2-RDM. The expression of these linear 
equations explicitly depends on which one-particle basis is employed. Thus, even though 
a trial 2-RDM is diagonal iV-representable, it is not necessary diagonal iV-representable 
on a different basis. Therefore we require its invariance, which can be included via usual 
semidefinite formulation. 

Our result yields a sufficient condition. Aside from significant theoretical importance, 
it is still impractical as incorporation of the diagonal iV-representability conditions in the 
RDM method has not yet been attempted, and the assumption that each subsystem being 
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the eigenstate of the number of electrons is not trivial as it is violated in practice [21]. 

This paper is organized as follows. In Section II, we briefly review the RDM method and 
7V-representability conditions. The main results are shown in Section III. The discussion 
and the conclusion are placed in Section IV. 

II. THE REDUCED-DENSITY-MATRIX METHOD AND THE N- 
REPRESENTABILITY CONDITIONS 

In this section, we briefly describe the RDM method and iV-representability conditions. 
Since the Hamiltonian involves one- and two-particle interactions, the Hamiltonian is thus 
expressed as follows: 

ij i\ilhh 

where a) and a are creation and annihilation operators, and v and w are one-particle and 
two-particle operators, respectively. The first- (1-RDM; 7) and second-order reduced density 
matrix (2-RDM; T) are respectively defined as follows: 

7} = hiatal*), 

and 

We can alternatively define using an ensemble state; J2 p w p \^/ p ){^/ p \, where w p is a non- 
negative value with unit sum, and { | ty p ) } is a complete orthonormal set of iV-particle space 
to define 1- and 2-RDMs. The ground state energy E g can be obtained using the 1- and 
2-RDMs as follows: 

E g = min(*|#|^) (1) 
= mmi^(^jat % |*) + - ^ {M^AA^h^) 

L ij 11123132 ) 

In this way, E g be calculated by 1- and 2-RDMs. We refer such an approach as "the 
RDM method". 



(2) 
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Until the 1960s, researchers tried to minimize the 2-RDMs using eqn. ([2]) instead of eqn. 
([1]). However, energy values obtained are far lower than the actual values, because their 
trial 2-RDMs are not derived from wavef unctions. The condition for 2-RDM to actually 
derived from a wavefunction is called the "iV-representability condition", which was first 



formulated by Coleman {2]. This condition states that if a given arbitrary 2-RDM satisfies 
this condition, and then it should have been derived from a wavefunction, but not otherwise. 
If it is non-iV-representable, no wavefunction is reduced to the 2-RDM. 

We denote Sn as the iV-representable set of the 1- and 2-RDM, and the RDM method is 
then expressed as follows: 



E 9 = ™ 7 {e^]+ £ 

\ ij tii2jija 



(3) 



Here are some of the trivial iV-representability conditions: 
fa) 1-RDM and 2-RDM are Hermitian, 



V = iVY T ili2 = (r^ 1 )* 

'j \H) 1 3132 \ tl»2 / ' 



(b) 2-RDM is antisymmetric, 



3132 ~ 3132 ~ 3231 ~ 3231' 



(c) the trace conditions are valid, or equivalently, the system is an eigenstate of the number 

of the electrons, 

i=l ij=l 

(d) a partial trace condition exists between 1- and 2-RDM, 

k=l 

From (c), we can prove (iV 2 ) — (iV) 2 = 0. In other words, the system is an eigenstate of 
the number of electrons. Of course, the above list is not exhaustive. For 1-RDM, the com- 
plete iV-representability conditions are quite easily obtained; a 1-RDM is ^-represent able 
if it satisfies the conditions stated above, and its eigenvalues are between zero and one 



12. 



tical 



28 L However, in the case of 2-RDM, obtaining a complete set of conditions is imprac- 
, |4[. Thus, finding a physically meaningful approximation or interpretation of the 



7V-representability conditions is crucial for the RDM methods. Some well-known necessary 
conditions are the P, Q ^j, and G defined by the positive semidefiniteness of P, Q, and 
G-matrix which in turn are defined as follows: 

G hn = (y\a{a l2 a} 2 a n \<U) h 0, 

where A >z denotes A is a positive semidefinite matrix; i.e. the eigenvalues of A are 
non-negative. For a system A ■ ■ ■ B, we denote £^"' B as 

&n B = {(r,7)|r satisfies the P, Q, and G conditions with trivial iV-rep. conditions from (a) to (d). } 
Using the P, Q, and G conditions, the RDM method can be formulated as a positive 

tin 

semidefinite programming problem pjj, [12|] . Then, the RDM method is represented as: 

In 2004, Zhao et al. [5|] implemented additional iV-representability conditions called the 
Tl and T2 conditions. Later, Braams et al. and Mazziotti [2^] derived T2' condition, 
which replaces T2 with slight enhancement. These conditions state that Tl, T2 and T2' 
matrices are positive semidefinite, and can also be written as a linear functional of 1- and 
2-RDMs. Therefore, we can still formulate as a standard type of semidefinite programming 
problem. In the RDM method, we usually obtain non-physical 2- RDM, in terms of having no 
wavefunctions that can be reduced to 2-RDM. Nevertheless, deviations seem to be negligible 
for the equilibrium geometry of molecules 0, Q, 14 - 18] . 



Other well-known iV-representability conditions include the diagonal iV-representability 



conditions 
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261 ] . These conditions state that if a trial 2-RDM V satisfies the diagonal iV- 
representability conditions, then there exist certain wavefunctions |^) which can be reduced 
to the diagonal part of original 2-RDM. 

rg = ^\ata}a jai \y) (5) 

Note that such a 2-RDM may have its ancestor wavefunctions |<3>) different from Ex- 
plicitly, we can construct |$) as follows: 

|$) =X}>Me*W, (6) 



where |<E>j) are the Slater determinants, <pi is an arbitra ry p hase factor and Wi are non 



negative numbers with unit sum (£\ Wi = l,Wi > 0) 25|. In this case, non-diagonal 
elements by |$) and j^) are likely to be different. 

Another important feature of the diagonal iV-representability conditions is that they can 
be written as a linear inequality of the diagonal elements of 2-RDM. Here we list some 
inequalities: 

7l>0 
l- 7 *>0 

r|>o 
7i - rg > o 
i- 7 i-7j + r|>o. 

Although all the inequalities can be written, expressing the complete set of diagonal repre- 
sentability conditions may be impractical, as the number of inequalities grow astronomically 



26l |. In this paper, we do not add any 



with size of the basis and the number of electrons 
new diagonal A^-representability condition, and we invesigate this physical nature. 

The expression of all these inequalities depends on which one-particle basis is employed. 
Thus, if a non A^-representable with a diagonal A^-representable 2-RDM is given, 7] > 
is satisfied for the specific choice of one-particle basis. In such a 2-RDM with a different 
one-particle basis representation, 7] > may be violated. Derivation of inequalities by 
the diagonal A^-representability conditions is independent of the choice of the one-particle 
basis. Therefore, we can generate more A^-representability conditions by applying a unitary 
rotation to the one-particle basis. In so doing, if the diagonal A^-representability conditions 



with unitary invariant are established [25[; e.g. if 7] > and 1— 7? > are unitary invariants, 
these two inequalities impose a constraint on eigenvalues of 7, which should be in [0, 1]. We 
can prove this as follows: since 7 is Hermitian, we can diagonalize 7. On the basis of 7] > 
7 is diagonalized, and the eigenvalues of 7 are greater or equal to zero. Next, when 1 — 7* > 
generates eigenvalues is smaller or equal to 1, the eigenvalue of 7 should be within [0, 1] on 
combining these two results. Interestingly, this is a complete A^-representability condition 
of 1-RDM. 

Apparently, every linear inequality has its unitary invariant form. If a given 2-RDM is 
diagonal A^-representable regardless of the one-particle basis employed, then we call such 



a 2-RDM as "unitary invariant diagonal A-representable" , and the conditions there of as 
"unitary invariant diagonal iV-representability conditions" . 



III. ON THE SIZE-CONSISTENCY 

In this section, after reviewing (i) the definitions of size-consistency and subsystem 1-, 
and 2-RDM and (ii) lemma and theorem for matrices, we obtain the main result to provide a 
proof (iii) using these lemma and theorem, if a variational set of 2-RDM satisfies the unitary 
invariant diagonal A-representability conditions and positive semidefinite conditions such 
as the P, Q, G, Tl and T2' . The 2-RDM of each subsystem derived is the eigenstate 
of the number of electrons, and the RDM method with this variational set is then size- 
consistent. Note that start formulating from using 1-RDM and cumulant of 2-RDM 341 ] . 



which are are additively separable quantities, do not solve the problem. Since considering 
A-representability condition on 2-RDM is same as using these quantities. As they contain 
exactly same information. 

A. The definition of size-consistency for the RDM method and subsystems 

Size-consistency may be defined as follows: assuming that A and B are two non- 
interacting systems. If a variational space is size-consistent, the total energy E(A ■ ■ ■ B) 
of the total system A - ■ ■ B is expressed as follows: 

E(A---B) = E(A)+E(B), (7) 



where E(A) and E(B) are energy values of the subsystems calculated separately [30(. 

Note that, as is already known, a variational method using 2-RDM with the P, Q and G 
condition as A-representability conditions is not size-consistent 2l|-|24]. 

First, let us consider the subsystems: 1- and 2-RDM of subsystems A (j A , T A ) and B 
(7 s and T B ) are defined using a set of one-particle basis allocated in each subsystem A or 
B from the whole system as follows: 

7 A * := 7*- where i,j correspond to the one-particle basis of A 

^ jih := ^hh where ix,i2, jx, 32 correspond to the one-particle basis of A 
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7 B ,- := 7* where i,j correspond to the one-particle basis of B 



j ■ 'J 

r B ''' 2 := Hf 2 where i\,i2, ji, J2 correspond to the one-particle basis of B. 

For other cases, 7 A * = is designated for either i or j corresponding to the one-particle 
basis of B, etc. 

We can define subsystem Hamiltonians H A and H B as well. First, one-particle operators 
and two-particle operators of subsystems A and B are expressed as follows: 

v A% j := Vj where i,j correspond to the one-particle basis of A 

w A *f 2 := w) 112 where «i, i 2? ji, H correspond to one-particle basis of A 

v B% j := V* where i,j correspond to the one-particle basis of B 

id 5 ' 1 ' 2 := w) lt2 where ii,i2,ji,j2 correspond to the one-particle basis of B. 

For other cases, w B jj 2 = is designated for convenience; if one of ii, i 2 , ji or j 2 corresponds 
to one-particle basis of A, etc. H A and H B are then express as follows: 



A J 1*2 t t 

w hh a h a i2 a h a n 



IJ IU2.71.J2 

and 



hh a ii a h a h a h ■ 



H B = J2v B ]a}a j+ £ 

ij hhjljl 

Note that definitions of T A ,'j A ,T B ,7 s , H A , H B are the same when we only consider the 
subsystem as the total system. 

Size-consistency within the RDM method requires: 



Ai\h T^Aiii2 

w 

h 



~..AlH2 -pAli 

hh hn I ' 



L ij hfajlj2 J 



and, 



E g {A---B)=E g {A) + E g {B). (8) 
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where J rA " B ,J rA and T B are sets of 1- and 2-RDMs satisfying certain iV-representability 
conditions for the total system A - ■ ■ B, for the subsystem A and B, respectively. 

Now, let us investigate the subsystems. To simplify the problem, and at the dissociation, 
each subsystem A and B should be the eigenstate of the number of electrons, thus we require 
each subsystem A and B to be the eigenstate of the number of electrons, i.e. Na and Nb 
electrons in the subsystem A and B, respectively, 

iy!=Afc, ^ A t= NAiN i~ 1} - 

i=l 



»J=1 



A variational space for the subsystem J^n a * s defined by T A ' B as follows: 



(r A , 7 A ) 



For T G F A - B , 

T iff := T llt ? where ii,i2,ji,j2 correspond to one-particle basis of A, 



r 



3 •' 



7 



N- 



I Sfc ^jk where i, j correspond to one-particle basis of A, 



andEi=i7^ = iV A , £ 



pAU _ N A (N A -1) 
i J=l *j 2 



In a similar manner, we can define J- B . Finally, we consider the size-consistency problem 
eqn. (jSJ) which is expressed as follows: 

min <W'y t - + V w hi2 T hi2 > = 

I ij ) 

^ Ai Ai , \^ Aiii 2 pAiii 2 

n«2iii2 

E*v;+ E w b tsj b ?s] 

ij i\ii3\h ) 



mm 



+ mm 



(9) 



Size-consistency disturbance happens when eq. (Q is violated. If only the necessary 
conditions are imposed; e.g. if we impose the P, Q, G conditions as iV-representability 



conditions, the left hand side decreases 



22|: 



E g {A---B)<E g (A) + E g {B), 

and if only sufficient conditions are imposed; e.g. if we use T A ' " B , Fn a {.F a '" b ) and 
jF b (F a " b ) as a set of 2-RDMs from the SDCI wavef unctions, the left hand side increases: 



E g (A---B)>E g (A) + E g (B). 
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B. The lemma and theorem for matrices 



In this subsection we show the lemma and the theorem of matrices. 

Lemma 1. If P is an nxn symmetric matrix and positive semidefinite, its mxm (n > m) 
principal submatrix A m is positive semidefinite. 



Proof 1. Proof is given in 



Theorem 2. Let A and B be n x n Hermitian matrices. If (WAU)u = (WBU)a, Vz G 
{1,2, ... ,n} is valid for all unitary matrices U , then A = B. 

Proof 2. Let us prove by induction on the dimension of the matrix. Ifn — 1, the result is 
trivial. Suppose it is valid for all matrices of size n — 1, consider all nxn unitary matrix U 
such that Ui n = U n i = 0(1 <i<n — 1). From the induction hypothesis, we conclude that 
Aij = Bij (1 < i < n — 1, 1 < j < n — 1), and A nn = B nn . In a similar manner, we can show 
that A = B except for the elements (l,n) and (n, 1), if Ui n = U n i = (2 < % < n). Now, 
let us call the first row of an arbitrary unitary matrix U by u\ and its last row by u n . Then, 
W(A-B)U = (A ln - B ln )u\u n + (A nl - B nl )ui Ul = (A ln - B ln )u\u n + {{A ln - B ln )u\u n )^ . 
In other words, from our hypothesis, the diagonal elements of this matrix should have zero 
real parts for any choice of unitary matrix U . Therefore, A\ n = B\ n = A' nl = B\ tl . ■ 



C. A set of approximate iV-representable 2-RDMs satisfying size-consistency 

We further investigate the nature of Fn a (£n B )- From Lemma 1, Y A G J^Na^n B ) ^ s 
positive semidefinite, since T is positive semidefinite and T A is the principal submatrix of 
T. Moreover, subsystem Q A and G A matrices defined as follows are also trivially positive 
semidefinite: 

Q A1 1a„ = Q % l l L where i\, 12,31,32 correspond to the one-particle basis of A, 

and 

G A *f. 2 = Gl 1 *? where i\, 12,31,32 correspond to the one-particle basis of A. 

The RDM method for the subsystem A then becomes: 



mm 



j* 3 3132 3132 



:F N a (. £ n" B )^ VA ^ A ) { ij hi 2 jlj2 
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Inconsistency arises since J~^ A (S§ "' B ) and are different, and is numerically shown to 



that Fn a (£n " b ) ^ Na |22j|- Let us investigate G condition in the subsystem A, 

hh ~ h< h Ji*2' 

where 5 is the kronker's delta. This G A matrix is apparently positive semidefinite. However, 
7^ carries some information from the total system A - ■ ■ B. If we use only the variables in 
subsystem A, we need to use 1-RDM 7^, taking the partial trace of Y A of subsystem A, not 
r of the total system A - ■ ■ B: 



k=l 



7 J jv,-1 



then we consider the G-condition of the subsystem A by defining G A matrix as follows: 

hh h ' h hh' 

Semidefiniteness of this G A should serve as the G-condition in subsystem A. However, 
7^ may not be consistent with j A . This inconsistency causes size- inconsistency in the 
RDM method. The RDM method finds a lower energy value by breaking condition (d), 
consistency between 1- and 2-RDM, of subsystem A (or B). Nevertheless, the following 
theorem establishes consistency between ■y A and j A . 

Theorem 3. // (possibly non-) N-representable 2-RDM T of the whole system satisfies 
unitary invariant diagonal N -representability conditions, and subsystem T A is the eigenstate 
of the number of electrons, then, 7 A * := N ^_ x Ylk=i ^ A fk = 1 A ) ^ s valid. 

Proof 3. It is enough to show that ^ A \ = ^ A \ for a specific choice on the one-particle basis. 
This is because the unitary invariant N -representability of V , and the expression of each 
linear inequality corresponding on the diagonal N -representability condition is independent 
of the choice of one-particle basis. Then, 7 A * = 7 A * is thus derived from Theorem 2. 

Let us prove 7 \ = r ) A \. From the diagonal N -representability condition, there exists |^) 
which reduces to the diagonal elements ofT, 



We can choose such |^) as 
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where |<& p ) are the Slater determinants, and w p is a non-negative value with unit sum 
CLj P \ w p\ 2 = -U while <f) p is an arbitrary phase factor I2aj . 

Since we employ separate one-particle basis in A and B, |$ p ) can be decomposed into 
|$^$^) := <S> \<&r) = \&p); where |$^) and ) are the Slater determinants in each 
subsystem. In this way, the total wave function \^f) can be written as: 



(10) 



where w pq is a non-negative value with unit sum (Yl pq \w pq \ 2 = 1), and <p pq is an arbi- 
trary phase factor. Now, let us show that an ensemble Yl pq w pq\^p) (^p I ^hat reproduces the 
diagonal elements ofY A . 



and 



ij 



-ma\a)a jai \^) 



pqrs 



pqrs 



'w pq e'™ y/w r ~e i(prs ($^$f | a^ajOi | $f ) 



'w pq e 



>w r 



-j2^p q (K\ a H a j a ^t) 



\><i 



tr ( -a\a)a jai ^w pq \<$>*)(^\ j = - ^ ^(^Ja^a^) 
\ pi J pq 



■y A \ can be derived from ^2 pq w pq \$ p ){$ A} 



1 A \ 



N A 



k \ p q / pq 



w pq (^\a\ ai \^). 
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Then, 7 \ = 7 \ can be shown as follows: 
2 



1 A \ 



N 

3 

2 



j 

1 E \y/^y/^** M ^** r '{^ \a\a)a jai \^ B q ) 



N 

jpqrs 



N 



1 \ w v^p® B M a ) a MK® B q ) 



2 



N- 
2 



N 



Vjr'eA.pg jeB,pq 



/ x ^ ^Aij _l N B ~ A j 
I i 



2 (yt^ + 

N - 1 ^ JJ 2 



AT — 1 V 



A physical interpretation in Theorem 3 requires each subsystem energy to be evaluated 
by the variables or information of each subsystem. 
Finally, we focus on Theorem 4. 

Theorem 4. The RDM method is size- consistent if the following conditions are satisfied; 
(i) the variational set of 2-RDM is unitary invariant diagonal N -representable, (ii) calcu- 
lated subsystems are the eigenstate of the number of electrons, (Hi) the 2-RDMs satisfy the 
arbitrary choice of positive semidefinite type N -representability conditions from conditions 
P, Q andG. 

Proof 4. Considering two non-interacting subsystems A and B being treated as a combined 
system A - ■ ■ B, the total energy of the combined system is now become the sum of the total 
energy of the subsystems; E(A ■ ■ ■ B) = E(A) + E(B). We define the approximate N- 

•A- 
N 



representable set of 1- and 2-RDMs of the total system £^"' B as follows: 



trivial conditions (a) to (d) with conditions P, Q, and G, 

and the unitary invariant diagonal N -representability conditions. 

14 



Then, we project the total system in subsystem A by: 



£ A 

C N A 



For E$~ B 

^ A 7ih = ^j\h w ^ ere h,i2, hi h corresponds to the one-particle basis of A, 
7 A * = jj where i,j corresponds to the one-particle basis of A, 
I^g = N A (N A - 1)12, and £ j7 ^ = N A . 



and the variational space can be defined by using only the variable in subsystem A or denoted 

=N A 

as £ A as follows: 



i A 

C N A 



T A il A 



trivial condition (a) to (d) with the P, Q, and G conditions, 
and the unitary invariant diagonal N -representability conditions. 



£ A . First, we show that £^ A C £ N . It is easy to 



It is sufficient to demonstrate that £^ A 
see that (a) Hermite condition, (b) anti-symmetricity and (c) trace conditions are satisfied 
for either 1- or 2-RDMs of £^ A - From Theorem 3, partial trace condition holds for the 
subsystems, thus validating (d) as well. Now let us construct the Q A and G A matrices as 
follows: 



Q A jih = Qfih w ^ ere hii2ijiih correspond to the one-particle basis of A, 



and 



G 



A tll2 



G™? where i\,%ii hi h correspond to the one-particle basis of A. 



Matrices P , Q A and G are all positive semidefinite, because from Lemma 1, all principal 
submatrices are positive semidefinite. Again, from Theorem 3, ^ A is derived from T A , thus 
matrices Q A and G A of subsystem A are expressed via T A and •y A . Furthermore, both sets 

— =A =A — 

satisfy the diagonal N -representability conditions to yield £^ A C £ Na - Next, £ N C £^ A , 
where we construct the 1 and 2- RDM 7 and V of the total system from subsystem 1- and 
2-RDM of A and B; -f A , j B , Y A and T B , respectively by: 



7 



1 A ) 

„Bi 



7"*- where i,j correspond to the one-particle basis of A, 



7"j- where i,j correspond to the one-particle basis of B, 
where otherwise. 
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and 




<Ahi2 
3132 



where %i, 12, 31,32 correspond to the one-particle basis of A 
where 11,12,31,32 correspond to the one-particle basis of B 
where i\,ji correspond to the one-particle basis of A, 
and 12,32 corresponds to the one-particle basis of B 
where 12,32 correspond to the one-particle basis of A, 
and ii,ji corresponds to the one-particle basis of A 
where 11,32 correspond to the one-particle basis of A, 
andi2,ji correspond to the one-particle basis of B 
where %2,3i correspond to the one-particle basis of A, 
andii,]2 correspond to the one-particle basis of B 




lryAh B%2 

2 ' 31 I 32 




2 ' 32 I 31 



r = < 




_1 Ail Bh 



2 I 32 1 31 



3132 



_l^,Ai2 Bh 
2 I 3i ' 32 



v 








otherwise. 



It is easy to show that T and 7 satisfy (a), (b), (c) and (d). The 2-RDMT has a structure 
containing three blocks; (i) where indices correspond to subsystem A, (ii) indices correspond 
to subsystem B, (Hi) indices correspond to both subsystems A and B. For condition P, it 
is sufficient to express these three matrices as positive semidefinite. Apparently, the first 
two blocks are positive semidefinite by their construction. The third block is always N- 
representable from the unitary invariant diagonal N -representability conditions (an explicit 
representation of the third block is derived from eq. / fTOj) and Theorem 2 employing an one- 
particle basis which diagonalizes subsystems A and B), thus T satisfies condition P. We can 
also show the positive semidefiniteness of matrices Q and G of the total system as well. Thus 

— =A — =A — 

(r,7) e £n" b , then £ N C S§ holds. Combining these two results, we obtain S Na = £^ A - 
Therefore, the total energy becomes the sum of each subsystem. ■ 

Here we have two corollaries. 

Corollary 5. Size- consistency is also satisfied by the RDM method employing 2-RDMs as 
long as the unitary invariant diagonal N -representability conditions with P , Q, G, Tl and 
T2' condition satisfied, assuming that each system is an eigenstate of the number of electrons 
at the minimum. 

Corollary 6. Size- consistency is satisfied by the RDM method employing 2-RDMs as long as 
the unitary invariant diagonal N -representability conditions with an arbitrary combination 



1(3 



of conditions P, Q, G, Tl, andT2' are satisfied, assuming that each system is the eigenstate 
of the number of electrons at the minimum. 

IV. DISCUSSION AND CONCLUSION 

We have showed size-consistency can be reestablished by the RDM method using uni- 
tary invariant diagonal iV-representability conditions in 2-RDM. This is an extension of the 
diagonal iV-representability conditions. 

We also assumed that in the RDM method, electrons are distributed to each subsystem 
appropriately This requirement looks stringent, as Van Aggelen et al. have showed that 
an incorrect distribution exists at the dissociation limit. However, we concur on the finding 
that the unitary invariant diagonal iV-representability conditions results in the correct dis- 
tribution of electrons at the non- interacting limit. Thus, we can omit this assumption. Note 
that we usually employ ensemble 7V-representability. In such a situation, even if we force to 
employ the complete ensemble to satisfy iV-represent ability, each subsystem would become 
a fractional number of electrons. For example, consider the dissociation limit of a diatomic 
isonuclear system with an odd number of electrons, e.g., Nj. It is not clear which N atom 
has seven electrons, and ends up as an ensemble of six to seven electrons. In this case, we 
need a pure iV-representability condition, or we need to apply certain perturbation to break 
the degeneracy. 

In practice, inclusion of the unitary invariant diagonal iV-representability conditions is 
not a simple task. A possible way is finding a positive semidefinite type iV-representability 
conditions corresponding to implying inequalities from the diagonal TV-representability con- 
ditions. For example, condition P implies Tfj > and it includes Tfj > for all one-particle 
basis representations. Apparently condition P is stronger than > since semidifinite- 
ness of P-matrix include unitary rotation of two-particle operators. Likewise, condition Q 
implies 1 — 7? — 7j + T^ > and condition G-condition 7* — Y 1 ^ > 0. Using anti-commutation 
relation of creation and annihilation operators 

a\a,j + aja\ = 5j, a^aj + OjOi = 0, and ajaj + ajaj = 0, 

where 5 is the Kronecker's delta, we can verify these relation by taking the diagonal elements 
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of Q and G-matrices. For Q-matrix as follows: 



gJJ = (tf|a il a ia at,a] 1 |tf) 

= (5 h 5 i2 - 5 il 5 i2 ) - (5 h Y 2 + cf-V- 1 ) + (<f.V- 2 + ^ 2 7 il ) - 2r il *? , 
v 31 n n n> K n <n n hi' y n hi 31 hi' 3132' 

and non-negativity of diagonal elements of Q-matrix implies following inequality, 

Q% = i-oj-7i-2rg>o. 

For G-matrix, we can show as follows: 

= («-2r-g)>o, 

and non-negativity of diagonal elements of G-matrix implies following inequality, 

G% = 7 « - 2rg > 0. 
Likewise, for Tl-matrix, we can show as follows: 

( T1 )}f/al = k% a \AA^ a h a h a h + a h a h a h a \AAi\^) 
= 2(S J2 T ilis -I- A^F* 2 * 1 -L- A-Jip^ia _i_ diyhh 

' IJ 3133 «3 Jail »1 3332 «3 3332 

_ri2pi2«i _ ^J3-p«3n _i_ xi3-p«3«2 _xii"p«3*i _ xh-phi2\ 
h 33j'i *2 jiia ' ii hh *a 3332 ii 33 ji/ 

+5 i3 5 jl 7 i2 - 5 i2 5 J1 7 i2 - i + 5 il ^ 3 7 i3 + 6 U 6 J - '•■ 

1 11 13 '32 II 13 '33 11 IJ '31 1 11 12 132 U IJ '31 



H 12 '33 1 «2 13 '33 12 13 hi »2 »3 '32 1 11 13 '31 



H 13 '32 1 <2 !3 '31 <2 13 '33 H «3 '31 1 H «3 '33 

S j3 5 J1 7* 1 - ^V- 3 + # 2 #y. 3 - <5 i2 cP<5 i3 - aWtf? 1 

12 *3 '32 11 IJ '32 11 12 'J3 II 12 13 11 12 13 

+5 j2 5 j3 5 jl + <P<Ptf? a + tf^tf?' 8 - S jl S j3 S j2 . 

ll 12 13 1 11 12 13 11 12 13 11 12 «3 

Then, the non- negativity of the diagonal part of Tl-matrix 

(Tl)f k = 2Tg + 21* + 21* - 7^ - 7l - + 1 > 
implies Condition VI: 1 - 7* - 7j - 7^ + 2rg + 2Tf k + 2Tf k . For T2-matrix, 

( T2 )hhh = (^l^^^^^^i + ^^^'i^ola^l*) 

= 2fT 3 r ili2 — A J1 r- ?3 * 2 -L ^2-pJ3*2 _|_ f3lp33il _ X32p33H\ 
33 3132 U «332 »1 1331 *2 1332 12 1331/ 

~ U l2 "ll "13 *2 H '13 ' 
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and the diagonal elements of T2-matrix is following: 

(T2)f k = 21* - 2rg - 2T* + j k k > 0, 

and its non-negativity implies Condition VII: 7* — 2r*j — 21^+21^ > 0, by just reordering the 
indices, of Weinhold- Wilson inequality, respectively. Non-negativity of T2' matrix implies 
Condition VI 



inequality 



and 7* > 26]. There exists such a correspondence for each condition and 
25| . By just adding all positive semidefinite type iV- representability conditions 
(e.g. P, Q, G, Tl and T2' etc.) is enough to satisfy the unitary invariant diagonal N- 
representability conditions. 



Interestingly, the unitary invariant diagonal iV-representability has already been proposed 



to be not only for a necessity but also for a sufficiency by Weinhold 29[. In fact, this has 



32j. However, as far as 



recently been proven at the 50th Sanibel symposium by Weiner 
the author's knowledge, a complete proof has never been in the literature. Of course if 
it is ture, size-consistency is automatically satisfied. The importance of our results is a 
direct proof using unitary invariant iV-represent ability conditions. Assuming the unitary 
invariant diagonal iV-represent ability conditions are necessary and sufficient, we can express 
the following corollary: 

Corollary 7. The complete set of positive semidefinite type N -representability conditions 
analogous to the diagonal representability conditions will be a complete set of the N- 
representability conditions. 

Such an extension is preferred because we need only the information from 2-RDM for 
the RDM method and not higher order conditions Q], we can still formulate the problem 
by semidefinite programming, and the approximate ^-represent able set equates asymptoti- 
cally exactly by adding iV-representability conditions from the diagonal iV-representability 
conditions. A positive semidefinite typ e condition seems to be much more stronger than the 
corresponding linear inequality jol. 33]. Thus we expect rapid convergence. 

Finding practical methods for incorporating the diagonal ^-representability conditions 
or the unitary invariant diagonal 7V-representability conditions is a challenge that warrants 
future investigation. 
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